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ABSTRACT 

The recent observations of type la supernovae suggest that the universe is accelerating 
now and decelerated in the recent past. This may be the evidence of the breakdown of 
the standard Friedmann equation. The Friedmann equation ~ p is modified to be 
a general form = g{p). Three model s with particu lar form of g(p) are considered in 
det ail. The supernova data pubhshed bv lTonrv et al] |2003) , D aly fc Diorgovskil l|200^ 
and lKnop et alJ ([2003^ are used to analyze the models. After the best fit parameters 
are obtained, we then find out the transition redshift zt when the universe switched 
from the deceleration phase to the acceleration phase. 

Key words: cosmological parameters-cosmology: theory-distance scale-supernovae: 
type la-radio galaxies: general 



1 INTRODUCTION 

The type la supernova observations suggest that the ex- 
pansion of the universe is speeding up rather than slowing 
down jPerlmuttt e r et al. 1998 ,1999; Garnavich pt al 199a 
iRiess et 81.11199^ : iTonrv et alJl2o63l: iKnop et alJl2003l) . The 
measurements of the anisotropy of the cosmic microwave 
background favor a fl at universe llde Bernardis et al J 120001 
i Hananv et all l200d : iBennett et alj l2003l : ISoergel et alJ 
I2OO3). The observation of type la supernova SN 1997ff at 
2; ~ 1.7 also supports that the universe wa s in the deceler- 
ation phase in the recent past llRiesjl200 J) . The transition 
from the deceleration phase to the acc eleration phase 
happene d around the redshift zt ~ 0.5 ijTurner fc RiessI 
[200a: .Dalv fc Diorgovskil l2003l) . A form of matter with 
negative pressure widely referred as dark energy is usually 
introduced to explain the accelerating expansion. The sim- 
plest form of dark energy is the cosmological constant with 
the equation of state pA = — PA. The cosmological constant 
model can be easily generalized to dynamical cosmological 
constant models, such as the dark energy model with the 
equation of state pq — luqpq, where the constant ujq 
satisfies —1 ^ ujq < —1/3. If we remove the null energy 
restriction ujq J5 —1 to allow supernegati ve ujq < — 1 , 
then we get t h e phantom energy models (lAlcanij 120031: 
ICaldwelll 120021: iKapUnghat fc Bridld l2003ll . More exotic 
equation of state is also possible, such as the Chaplygin 
gas model with the equation of state p — —A/p and the 
generalized Ch aplygin gas model with the equation o f state 
p = -A/p" llKamenshchik. Moschefla fc PasauieJ l200lt 
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also called the quintessence field. The energy density of the 
quintessence field must remain very small compared to that 
of radiation and matter at early epoches and evolves in a 
way that it started to dominate the universe around the 
redshift 0.5. The re are other forms of dark energy, like the 
tachy on fi eld l|A rmendariz-Picon. Dam our fc Mukhanovl 
\ i99d. IPadmanabhan fc Choudhurvl 120012 003: 
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Although dark energy models are consistent with cur- 
rent observations, the nature of dark energy is still a 
mystery. Therefore it is possible that the observations 
show a sign of the breakdown of the standard cosmology. 
Some alternative models to dark energy models were pro- 
posed along this line of reaso ning. These models are mo- 
tivated b y brane cos mology l*Dvali. Gab ada(j| zg_fc Po rratil 
[2003; Dva h fc Turiieil 12003: Chung fc Freesd 119991) . In this 
scenario, our universe is a three-brane embedded in a five 
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boundary 
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dimensional spacetime. The five dimensional action is 

-1 
2^ 

where TZ. is the Ricci scalar in five dimensions, G is the 
five dimensional metric determinant, K5 is the five di- 
mensional Newton's constant. Sorb is the orbifold action, 
^boundary represents the boundary action and Sqh is the 
Gibbons-Hawking boundary terms. In these models, the 
usual Friedmann equation = SnGp/S is modified to 
a general form = g{p) a nd the universe is composed 
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is the brane cosmology with g{p) ^ p + p . In order to re- 
tain the success of the standard cosmology at early times, 
we require that the modified cosmology recovers the stan- 
dard cosmology at early times. So g{p) must satisfy g{p) ~ p 
when p ^ po, where po is the current matter energy density. 

For a spatially fiat, isotropic and homogeneous universe 
with both an ordinary pressureless dust matter and a min- 
imally coupled scalar field Q sources, the Friedmann equa- 
tions are 

2 _ 8tvG , 
3 

(Pm + PQ + 3pQ), 



H 



-(Pm +Pq), 



(1) 



a 3 

PQ + 3-H'(pQ +pq) = 0, 



(2) 
(3) 

where dot means derivative with respect to time, pm — 
Pmo(fflo/cf)'^ is the matter energy density, a subscript means 
the value of the variable at present time, pq = /2 + V{Q), 
PQ = q2/2 - V{Q) and V{Q) is the potential of the 
quintessence field. The modified Friedmann equations for 
a spatially fiat universe are 

(4) 
(5) 
(6) 



= H^g{x), 



- = Hog{x) —g (x) 

a Z \ p 



p + 3H{p + p) = 0, 

where x = SttGp/SHq = xo{l + z)^ during the matter dom- 
inated epoch, 1 + z — ao/a is the redshift parameter, g{x) = 
X + ■■■ is a general function of x and g'{x) — dg{x)/dx. 
Note that the universe did not start to accelerate when the 
other nonlinear terms in g{x) started to dominate. To re- 
cover the standard cosmology at early times, we require that 
g{x) « X when x ^ xq. For the matter dominated fiat uni- 
verse, p — Pm and p — Pm = 0- Let QmO = SttGpo/SHq, 
then xo = fimo, g{xo) = 1 and x — f2mo(l + 2)^ during the 
matter dominated era. 

The luminosity distance di, is defined as 



di,{z) = ao(l + 2) 



dt' 



Ho 



g ^^^[Qmoii- + uy^]du. 



(7) 



The apparent magnitude redshift relation is 

m{z) = M + 51ogi(, dL(2)-h25 = A1+51ogu)©L(2) 



AH-51ogio <^ (1 + 2) / g-'-'^[nmo{l + uf]du\ 



(8) 



where I'l(2) = Hodi^(z) is the "Hubble-constant-free" lu- 
minosity distance, M is the absolute peak magnitude and 
M = M — Slogj^Q Ho -\- 25. M can be determined from the 
low redshift limit at where T>i,{z) — 2. The parameters in 
our model are determined by minimizing 



E 



[rriobsizi) — in{zi)f 



(9) 



where ai is the total uncertainty in the observations. The x - 
minimization procedure is based on MINUIT code. We use 
the 54 supernova data with both the stretch correction and 
the host-gala xy extinction corre ction, i.e., the fit 3 super- 
nova data by iKnop et al .' ^20oS) , the 20 r adio galaxy and 
78 supernova data bv IPalv fc Dior gqvskil (i200^, and the 
supernova data by iTonrv et^a l. ( 2003') to find the best fit 
parameters. In the fit, the range of parameter space for M 
is = [—3.9, 3.2], the range of parameter space for firno is 
QmO = [0, 4]. 

The transition from deceleration to acceleration hap- 
pens when the deceleration parameter q — —ci/aH^ — 0. 
From equations 101 and JSJl, we have 

5[!^nr0(l+2T)^] = ^ a^rO ( 1 -f^ 2t) V [f^mO (1 2t )^] , (10) 

3 

qO = -fimOp'(fimo) — 1. (11) 

To compare the modified model with the dark energy model, 
we make the following identification 



9{.x) - 



(12) 



2 CHAPLYGIN GAS MODEL 

The chaplygin gas model p = —A/p" in the framework of 
alternative model to dark energy is 

g{x) = X + ncio[As + As)[x/Q.mof]^''', 

where SIqo = 1 — fimo, P = 1 -I- a and As = 
(87rG/3i?of2Qo)'^ A. The a = 1 model is motivated by a d- 
brane in ci-|-2 spacetime. Since g'{x) = 1-\'Q.qo(^- — A s)[A a -\- 

(1 - As){x/nmofY"'-^{x/nmof, SO 



20, 



QO 



(1 + Z,^of[As + (1 - ^.)(1 + Z.^of^]'-'/" 



= A,--(l-A)(l + 2,=o)"', (13) 



1 3 



(14) 



go < gives that As > [1 — f2mo)~^/3. To retain the success 
of the standard model at early epoches, we require g{x) ~ x 
when X ^ 1. In other words, we require As ~ 1. Therefore, 
we have the following constraints 

< As <1, (15) 
As ~ 1. (16) 



The best fits to the 54 supernovae by iKnop et al.l 

are SlrnO = [0, 0.44] centered at almost zero. As — [0.99, 1] 
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centered at almost one and /3 — [1, 32.3] centered at 22.0 
with ~ 43.9. The best fi ts to the 98 radio galaxy and 
supernova data compiled bv lPalv fc Diorgovskil l)2003fl are 
n^o = [0, 0.27] centered at 0.26, As = [0.65, 1] centered 
at 0.97 and /3 = 1 with = 87.8. The best fits to the 
172 supernovae wi th redshift z > 0.01 and Av < 0.5 mag 
jTonrv et al.ll2003l) are fl^o = [0, 0.3] centered at 0.1, = 
[0.99, 1] centered at almost one and (3 = [1.5, 23.7] centered 
at 16.1 with — 169. 5. Th e best fits to t he 194 supernovae 
by iTonrv et al.1 ll2003l) and iBarris et all (Eooi) are Q^o = 
[0, 0.36] centered at 0.19, As = [0.99, 1] centered at almost 
one and p = [1.3, 25.2] centered at 14.9 with = 195.1. 
The best fits to all the data combined are fimo ~ [0, 0.34] 
centered at 0.21, As = [0.99, 1] centered at almost one and 
(3 = [1.0, 21.2] centered at 13.7 with x^ = 331.3. From the 
above results, we conclude that the generalized Chaplygin 
gas model tends to be the A model. 



3 GENERALIZED CARDASSIAN MODEL 

The model is 

<?(a;) = s[l + fe°<"''^]'''", 

where B = - l)/r2°'o""'\ q > and n < 1 - 1/3(1 - 

n^o). When n = 0, p(x ) = B^ '^d + which is the 

case studied by iFreesd ll2003l) . For the special case a = 1 
and n = 0, g{x) = x + B which is the standard cosmology 
with a cosmological constant. From a purely phenomenolog- 
ical point of view we may think that gravity is modified in 
such a way that acceleration kicks in when the energy den- 
sity approaches a certain value ll|Carro ll 2003). The model 
is motivated from a three-brane located at the Z2 symme- 
try fixed plane of a five dimensional spacetime. Chung and 
Freese showed that if one parametrizes the Hubble rate in 
terms of the brane energy density, then Cardassian model is 
derived with suitab le choice of the five d imensional energy 
momentum tensor llChung fc Freeselll999^ . The generalized 
Cardassian model gives 

g'{x) = [1 + Ba:"<"-'']'/" 

+ {n - l)Bx"'-"-'''>[l + Bx"'-"-^^]^/"'\ (17) 

Combining equation 11711 with equations (IIOII and lllll . we 
get 

l + ^T = [(t^-^-l)(2-3n)]i/^"(i-"', (18) 

go = i + |(n-i)(i-ns;o). (19) 

If we think the generalized Cardassian model as ordinary 
Freidmann universe composed of matter and dark energy, 
we can identify the following relationship for the parameters 
in the Cardassian and quintessence models 

t^QO = : ^ • 

J- — i'mO 

There are four parameters in the fits: the mass den- 
sity QrnO, the parameters n and a, as well as the nuisance 
parameter A4. The range of parameter space explored is: 
n = [-10, 0.6 6] and a = (0, lO"^]. The best fits to the su- 
pernova data JPalv fc Diorgovskil l2003l : iTonrv et alJ 120031: 
iKnoD et alJ l2003^ generally give very large a > 100, so 
B « ^^mo" and the transition redshift is weakly dependent 



n and n Confidence Contour with a=1 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



Figure 1. The 70%, 95% and 99% confidence contours of f2mO 
and n in Cardassian model from the 54 supernovae given by Knop 
et al. (2003) 

on a. Furthermore, x^ changes very little when a changes 
over a fairly large range. In other words, the generalized Car- 
dassian model differs little from the Cardassian model with 
a — 1. So we will discuss the Cardassian model in more 
detail below. 



3.1 Cardassian Model 

The Cardassian model is the special case a = 1 of the 
generalized Cardassian model. This model is equivalent to 
the dark energy model with a constant equation of state 
PQ = ^QpQ in the sense of dynamical evolution. The equiv- 
alence is provided by n = 1 -I- uqo and the equivalent dark 
energy potential is V{Q) = j4[sinhfc(Q/Q: -I- C)]~" with 
Q = —2 — 2/( n — 1) . The best fits to the 54 supernovae 
bv lKnop et aO (120031) are 0.^0 = 0.56to°il, n = S-Gtli 
and x^ ~ 43.73. The flmo and n contour plot is shown in 
figure Q 

The be st fits to the 98 radio galaxy and supernova data 
compiled bv lDalv fc Diorgovsk i ( 2003.) are flmo = O.Uto'H, 
n — 0.261q gj^ and x^ = 87.45. The contour plot is shown in 
figure 121 

The best fits to the 1 72 supernovae with redshift z > 
0.01 and Ay < 0.5 mag jTonrv et alJ |2003^ are r2mO = 
0.48t;;;?^, n = -l.2t\ i and = 171.4 Th e best fits to 
the 19 4 supernovae bv lTonrv et"ai] ^2003^ and lBarris et alJ 
J2004D are fimo = 0.5llo°6, n = -1.2+i;g and x^ = 196.7. 
The contour plot is shown is figure El T he plot agrees well 
with t he fig ure 13 in lTonrv et al.l i2003r and the figure 6 in 
iFrithI ll2004t) . 

The above results are summarized in tabled Combining 
the above results, we find that flmo ~ [0, 0.61] centered at 
0.45 and n = [—3.1, 0.5] centered at —0.6 at 99% confidence 
level. The 99% contour plot is shown in figure^] Take Qmo ~ 
0.3, we get zt = 0.35, go = —3.1 and u)qo = —3.44 when 
n = -2.44; zt = 0.57, go = -0.24 and loqo = -0.7 when 
n = 0.3. These results are consistent with those obtained 
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Table 1. Best fits to Cardassian model 




Fit 




mO 








# 


70% 


99% 


70% 


99% 




1 




'^■00_o.56 


■^■"-3.4 




43.73 


2 


0.14^°;?^ 


14+°-48 

'^•^^-0.14 




0.26«:- 


87.45 


3 


0.48i°;f8 




-1 2+^-1 


-1 2+l'^ 
-^•^-4.7 


171.4 


4 






-1 2+^-^ 


-i-2ll:«9 


196.7 


5 


Q 45+0.10 


45+0.16 
'^■^^-0.45 






332.0 



Fit 1 is tile fit to the 54 supernova data from Kno p et al. |20(^ , 
fit 2 is the fit to the 98 data points from IPalv fc Dj orpjovski 
J2003l ). fit 3 is the fit to the 172 supernova data fromlTon^^^al ' 



1 2003), fit 4 is the fit to the 194 supernova data from lTonrv et al 
12003) and Bar ris et al.. 12004*) and fit 5 is the fit to the above 
data combined. 



Figure 2. The 70% and 99% confidence contours of f2niO and n 
in Cardassian model from 20 radio galaxies and 78 supernovae 
compiled by Daly & Djorgovski (2003) 
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^1 99% (Tonry et al) 
2 Z 99°''° PtJiy el al) 
^ ; ;: 99% (Knop et al) 




Figure 4. The 99% confidence contours of CmO and n in Cardas- 
sian model from the sample data in Tonry et al. (2003), Daly & 
Djorgovski (2003) and Knop et al. (2003) 



Figure 3. The 70%, 95% and 99% confidence contours of f2m0 
and n in Cardassian model from the 172 supernovae with z > 0.01 
and < 0.5 mag listed in Tonry et al. (2003) 

by Zhufc Fujimoto (2003a.b.c'). IZhu. Fuiimoto fc Hd ll2003l) 
and lSen fc Sen. (.2003^ . 



4 MODEL 3 

The last mod el we would like to consider is pix) = 
[a + Va'^+xf JPeffavet et alJl2002l: iDvali fc Turneill2003l) . 
where a = (1 — Om o)/2. This model arises from the 
brane world theory by iDvali. Gabada dzc fc Pprrat i (2003) 
in which gravity appears four dimensional at short distances 
while modified at large distances. For this model, we find 
that the equivalent dark energy equation of state parame- 
ter ujQo, go and the transition redshift zt from decelerated 



expansion to accelerated expansion are 

t^QO 

qo 



-1/(1 + a^o), 

2f2mO — 1 



mO 



l + ZT 



2(1 - Omo) 

f2mO 



1/3 



(20) 
(21) 

(22) 



Applying the 54 supernova data with host-galaxy extinction 
correction (Knop_ct al. 2003), we find that r^mo ~ 0.19+q q5 
and = 45 .71. The 20 radio galax y and the 78 su- 
pernova data ijPalv fc Diorgovskil l2003l) give the best fit 
fimo = 0.18 ± 0.03 and = 87.6. The best fit from the 
172 supernovae with redshift z > 0.01 and < 0.5 mag 
is n^no = 0.nt.o°t a.nd = 175.2. I f we use the 194 
supern ovae given by iTonrv et alJ ll2003l) and iBarris et alJ 
J200'tf) . we find that n^o = 0.22l°;°* and x^ = 200.6. 
The above results are summarized in table |5| Combining 
the above results, we get QmO = 0.20 ± 0.02 at la level or 
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Table 2. Best fits to model 3 



Data # n,-no 2t ^^qo 

source Icr 3(T 

Knop 54 0.19loo5 0-19to.i2 45.71 0.90 -0.84 

Daly 98 0.18 ±0.03 O.lS+g g^ 87.6 0.95 -0.85 

Tonry 172 0.17+°;„3 0.17+oJg 175.2 1.0 -0.85 

Barris 194 0.22+0^^ 0.22l:oJg 200.6 0.77 -0.82 

Combined 346 0.20 ± 0.02 0.20^0.06 334.6 0.86 -0.83 
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f2mO = 0.201 Q at 3a level. If we take OmO ~ 0.15, then we 
have qo — —0.61 and zt = 1.13. If we take flmo ~ 0.21, then 
we have go = —0.48 and zt ~ 0.81. If we take QmO ~ 0.28, 
then we have qo — —0.34 and zt = 0.55. Thes e results are 
consistent with those obtained bv lDeffavet et al. 12002.) . 



5 DISCUSSIONS AND CONCLUSIONS 

A general function g{x) of the ordinary matter density was 
used to explain the current accelerating expansion of the 
universe. In this model, no exotic matter form is needed. 
This approach is equivalent to dark energy model building 
approach in the sense of dynamical evolution of the uni- 
verse because we can map the modified part of energy den- 
sity to dark energy. The function g{x) satisfies the following 
conditions: (1) g{xo) — 1; (2) g{x) ^ x when z 3> 1; (3) 
g{xo) > 3xog'{xo)/2, where x = r2mo(l + z)'"^ + f2ri)(l + z)". 
The generalized Chaplygin gas model tends to take As — 1 
which becomes the dark energy model with a cosmological 
constant. Therefore the generalized Chaplygin gas model is 
disfavored in the framework of alternative models although 
it is a viable dark energy model. Unlike the gravitational 
lensing constraint l)Dev. Alcaniz fc Jainl l20o3l. the super- 
nova data do not provide tight constraint on the generalized 
Cardassian model. A fairly large range of parameters from 
the generalized Cardassian model are consistent with the 
supernova data. For the Cardassian model, the supernova 
data give QmO = [0, 0.62] and n = [-3.11, 0.3] at the 99% 
confidence level. If we have better constraint on the tran- 
sition redshift zt, then we will be able to distinguish the 
Cardassian model from the A model because the A model 
is the special case n = 0. For the model 3, we find that 
Qmo = [0.12, 0.28] with 99.7% confidence. Only the upper 
limit gives zt ~ 0.5. 
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